Casimir effect in a superconducting cavity and the thermal controversy 
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One of the most important and still unresolved problems in the field of dispersion forces, is that of 
determining the influence of temperature on the Casimir force between two metallic plates. While 
alternative theoretical approaches lead to contradictory predictions for the magnitude of the effect, 
no experiment has yet detected the thermal correction to the Casimir force. In this paper we show 
that a superconducting cavity provides a new opportunity to investigate the problem of the thermal 
dependence of the Casimir force in real materials, by looking at the change of the Casimir force 
determined by a small change of temperature. The actual feasibility of the proposed scheme is 
briefly discussed. 
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I. INTRODUCTION 



Dispersion forces have been the subject of intense the- 
oretical and experimental investigations for a long time. 
This is not surprising, because these weak intermolecu- 
lar forces have a truly pervasive impact, from biology to 
chemistry, from physics to engineering lj . It may there- 
fore come as a surprise to know that there still exist, in 
this well established field, unresolved problems of a fun- 
damental character. This is indeed the case with respect 
to the problem of determining the van der Waals-Casimir 
interaction between two metallic bodies at finite temper- 
ature [3, y, Q. As of now, people simply don't know how 
to compute it, and the numerous recent literature on this 
subject provides contradictory recipes, which give widely 
different predictions for its magnitude @, S S|- Apart 
from its intrinsic interest as a problem in the theory of 
dispersion forces, addressing this problem is important 
because many experiments on non-ncwtonian forces at 
the submicron scale use metallic surfaces at room tem- 
perature [5j. We remark that, while the Casimir pres- 
sure has now been measured accurately in a number of 
experiments [a, Q, there is at present only one experi- 
ment that detected the influence of temperature on the 
Casimir-Polder interaction of a Bose-Einstein condensate 
with a dielectric substrate Q- There are however no ex- 
periments so far that detected the temperature correction 
to the Casimir force between two metallic bodies, and in- 
deed there are at present several ongoing and planned ex- 
periments to measure it [8| . Since the theoretical solution 
of this difficult problem is still far, it is of the greatest im- 
portance to do experiments to probe the effect. A recent 
accurate experiment using a micromechanical oscillator 
Q seems to favor one of the theoretical approaches that 
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have been proposed, but this claim is disputed by other 
researchers [3]. 

In this paper, we propose a new experiment using su- 
perconducting Casimir devices, that could help to find 
an answer to the question. Superconducting cavities are 
valuable tools to explore important aspects of Casimir 
physics. Some time ago we proposed to use rigid su- 
perconducting cavities to obtain the first direct measure- 
ment of the change of Casimir energy that accompanies 
the superconducting transition [9(. Here, we argue that 
it might be possible to shed some light on the contro- 
versial problem of the thermal Casimir effect in metallic 
systems, by measuring how the magnitude of the Casimir 
force changes as an effect of a change of temperature of 
the superconducting apparatus. Indeed we show that 
certain extrapolations to the superconducting state of 
existing alternative approaches to the thermal Casimir 
effect for normal metals, that have appeared in the re- 
cent literature, lead to strikingly different predictions for 
the magnitude of the effect. A distinctive feature of the 
proposed scheme is that, being a difference force measure- 
ment, it might achieve a better sensitivity than absolute 
measurements of the Casimir force. The idea of exploit- 
ing difference force measurements to probe the features 
of the thermal Casimir effect is not new indeed. Already 
in Ref.[10] the possibility of observing the difference in 
the thermal Casimir force in a normal (i.e. non supercon- 
ducting) metallic system, at two different temperatures, 
was considered, opening up a new opportunity for the 
investigation of the thermal Casimir effect in the submi- 
cron separation range. 

The plan of the paper is a follows: in Sec. II we briefly 
review the thermal Casimir effect, and explain in some 
detail the theoretical problems that arise in the case of 
metallic cavities. In Sec. Ill we describe our supercon- 
ducting Casimir devices and explain the general ideas 
on which our proposed experiment is based. In Sec IV 
we discuss the important issue of the possible prescrip- 



tions for the reflection coefficient of quasi-static magnetic 
fields by a superconducting plate, and in Sec V we apply 
these prescriptions to the computation of the predicted 
change in the Casimir force, in a plane-parallel system, as 
the temperature of the superconductor is varied. In Sec 
VI we discuss the experimentally important plate-sphere 
geometry, and finally in Sec. VII, we present our conclu- 
sions and briefly discuss the prospects of the proposed 
scheme of being experimentally feasible. 



II. THE THERMAL CASIMIR EFFECT IN 
METALLIC CAVITIES 

All current approaches to the thermal Casimir effect 
are based on the theory of the Casimir effect developed 
long ago by Lifshitz [11[ > on the basis of Rytov's general 
theory of e.m. fluctuations [12 |. According to this the- 
ory, the Casimir pressure P between two plane-parallel 
(possibly stratified) slabs at temperature T, separated by 
an empty gap of width a, is: 
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where the plus sign corresponds to an attraction between 
the plates. In this Equation, the prime over the Z-sum 
means that the I — term has to taken with a weight 
one half, T is the temperature, k^ denotes the projec- 
tion of the wave-vector onto the plane of the plates and 
1i = V k ± +£f/ c2 ' where & = 2irk B Tl/h are the Mat- 
subara frequencies. The quantities r„ = r« (i£;,kj_), 
with j = 1,2, denote the reflection coefficients for the 
two slabs, for a-polarization (for simplicity, we do not 
consider the possibility of a non-diagonal reflection ma- 
trix) , evaluated at imaginary frequencies i£i . It is of the 
outmost importance at this point to make a remark on 
the range of validity of Eq. |T]): Lifshitz original deriva- 
tion assumed that spatial dispersion is absent, so that 
the slabs can be modeled by permittivities e^(u>) and 
/i^'(oj) depending only on the frequency u>, and the re- 
flection coefficients r^ («£;, kj_) have the familiar Fresnel 
expression. Nowadays it is well known [13| that Eq.fJJ 
is actually valid also for spatially dispersive media (sep- 
arated by an empty gap), like metals at low temperature 
and obviously superconductors, provided only that the 
appropriate reflection coefficients are used. 

We now come to the main topic of the present paper, 
i.e. the application of Lifshitz formula to non-magnetic 
metals (we set fJ-=l from now on), which is the experi- 
mentally important case. We suppose for simplicity that 
the plates are made of the same metal. For our pur- 
poses, it is convenient first to write Eq. |T]) as a sum 
of three terms, P , P and Pi, which include, re- 
spectively, the TE zero mode, the TM zero mode and all 
non- vanishing Matsubara modes: 



We consider first Pi (a, T) . This term poses no particular 
problems. At room temperature (300 K), the frequency 
£i of the first non- vanishing Matsubara mode that enters 
in Eq. (fT]) is around £i = 0.159 eV/Ti, and belongs to 
the IR region of the spectrum. On the other hand, the 
presence of the exponential factor exp(2 aqi) cuts off the 
Matsubara modes with frequencies greater than a few 
times the characteristic cavity (angular) frequency lo c = 
c/(2a). For the relevant experimental separations in the 
range from 100 nm to one or two microns, the frequency 
to c falls somewhere in the range from the near UV to 
the IR. In the frequency range extending from £i to a 
few times uj c , spatial dispersion is negligible and one can 
use the Fresnel formula for the reflection coefficients, in 
terms of the permittivity e(i£„). The values of e(i£, n ) can 
be obtained, using dispersion theory, from optical data 
of the plates (see however the footnote before Eq. (JTDJ) 
below) , or alternatively can be computed using analytical 
models, like for example the plasma model of IR optics: 
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where Qp is the plasma frequency. If needed, one can 
possibly augment the above simple model by oscillator 
contributions accounting for interband transitions that 
may become sizable at the higher frequencies j5j- If de- 
sired, one can also take account of dissipation, via the 
inclusion of a suitable set of relaxation frequencies in the 
various terms. For the experimentally relevant separa- 
tions a, the inclusion of all these corrections changes a 
bit the magnitude of Pi (a, T) , and poses no problems at 
all. 

We consider now the two terms in Eq. @ that have 

I = 0, namely P (TE) (a,T) and P (TM) (a,T) (we shall re- 
fer to these terms as zero modes). Obviously, a van- 
ishing complex imaginary frequency is the same as a 
real vanishing frequency, and therefore the zero modes 
involve the reflection coefficients of the slabs for static 
cm. fields. Let us examine first P (a,T). A zero- 
frequency TM mode represents a static electric field. 
Since a metal expels an electrostatic field, we clearly have 
ttm(0, kj_) = 1. Inserting this value of ttm into the I = 
term of Eq. (fT]) we see that in metals P (T) attains 
its maximum possible value 
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P(a,T) = P 



(TE) 



(a,T) + P^™\a,T) + P 1 (a,T). (2) 



where £(3) ~ 1.20 is Riemann zeta function. 

Consider now the TE zero mode P (a,T). Sur- 
prisingly, this single term has become the object of a 
controversy, in the last ten years, and no consensus has 
been reached yet about its correct value. The whole crisis 
of the thermal Casimir effect for metallic bodies in fact 
stems from this very term, and it is therefore useful to 
review briefly what the problem is. A zero-frequency TE 
mode represents a static magnetic field. Since non mag- 
netic metals are transparent to magnetic fields, it seems 



natural to suppose that we should take 
r TE (0,k i )| Dr =0. 



(5) 



When this value is inserted in Eq. ([T]) one finds easily 
that the TE zero mode gives a vanishing contribution to 
the Casimir pressure: 



P 
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This result is dubbed in the Casimir community as the 
Drude- model prediction [14 |. because rxE(0,k^) = is 
the value that is obtained by inserting into the Fres- 
nel formula for the TE reflection coefficient the familiar 
Drude extrapolation of the permittivity of a metal to low 
frequencies: 



e D (uo) = 1 



Q, 2 p/[uj(u + i-y)] . 



(7) 



In this formula, 7 is the relaxation frequency accounting 
for ohmic conductivity. Despite its physical plausibility, 
this Drude approach has been much criticized in the re- 
cent literature. The first problem that was encountered 
is theoretical, as it was shown that, in the case of per- 
fect lattices, the Drude model value of the TE zero mode 
is apparently inconsistent with Nernst theorem [lj| . The 
problem is subtle though, and the actual existence of this 
inconsistency is still disputed. For a thorough discussion 
of different points of view on this issue we refer the reader 
to the recent review Q (see also Ref-jlq]). In addition, 
and perhaps more importantly, it has been claimed that 
the Drude model is also inconsistent with recent experi- 
ments Q . The solution that has been proposed Q to cure 
both difficulties is very puzzling: in evaluating the low- 
frequency contributions of Eq. |[T]). and in particular the 
TE zero mode, instead of the physically plausible Drude 
model, one should use the plasma model of IR optics, 
Eq. ([3]) , extrapolated without modifications all the way to 
zero frequency. For the TE zero mode, this prescription 
leads to the following expression for the reflection coeffi- 
cient: 



7- T E(0,ki)| p i 
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and we shall refer to this as the plasma model prescrip- 
tion. For the relevant values of k± ~ l/(2a), and for typ- 
ical values of fip ~ 9eV/h~ (corresponding to gold), it is 
easy to verify that tte(0, kj_)| p i is close to one. For exam- 
ple, for a = 200 nm, we obtain ?"te(0, l/(2a))| p i = 0.90. 
Since tte(0, k_i_)|pi 7^ 0, Pq '(a, T) does not vanish any- 
more, and we let P (a,T; dp) | p i its magnitude, as de- 
termined by the plasma prescription. For lu c /Qp <C 1, 
P (a, T;fip) | p i has the following expansion: 



P (TE) (a,T;tt P )| pl 
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where <5 = c/Vtp. We see that, for sufficiently large sep- 
arations a, the magnitude of P (a, T\Vtp) \ v \ is com- 
parable to that P (a, T). Recalling that the Drude 

prescription predicts instead P |di — 0, we see that 
the predicted Casimir pressures for the two prescriptions 
differ by a quantity AP(a,T) = P(a,T)\ Dl - P(a,T)\ pl 
HI: 



(TE), 



AP(a,T) = -P ^(a,T;fi F )| p i. 



(10) 



(9) 



The minus sign in the r.h.s. means that the Drude model 
predicts a smaller Casimir pressure than the plasma 
model. 

The disagreement between the two prescriptions 
is most evident at large separations. Indeed, for 
kBTa/(hc) ^> I, Pi (a, T) becomes negligible and the en- 
tire Casimir pressure is given by the zero modes contribu- 
tion only @,[|. Since in this limit P (TE) (a,T; fi P )| pl ~ 

P (T) , we see at once that the pressure predicted by 
the Drude prescription is one-half that predicted by the 
plasma model prescription! Obviously then the straight- 
forward way to clarify the problem would be to measure 
the Casimir pressure at large separations, say around 5 
/im, for room temperature. Unfortunately, at such large 
separations the Casimir pressure is very small, and as a 
result all attempts made so far to measure it have been 
unsuccessful. 

At the time of this writing, the most significant exper- 
iments are the ones reported in [5| , in which the Casimir 
pressure was measured for separations in the range from 
160 to 750 nm. For such small separations, the quan- 
tity AP(a, T) is only a small fraction of the absolute 
Casimir pressure (for gold plates at room temperature, 
about two-three percent at separations of around 200- 
300 nm, where the experiment in the second of Refs.Q 
was most sensitive). Observing such a small correction 
via an absolute Casimir measurement is very difficult, 
because one needs to measure with high precision the 
plates separation, and control a number of systematic er- 
rors, like residual electrostatic attractions [17} , roughness 
of the plates etc. After an accurate analysis of these pos- 
sible sources of systematic errors, the authors of Ref.Q 
conclude that the Drude approach is inconsistent with 
the experiment at high confidence level. However, this 
statement has been disputed by other researchers (see the 
second of Refs.[2]). In this regard, it has been remarked 
recently [I8J that obtaining an accurate theoretical pre- 
diction, at the percent level, for the magnitude of Casimir 
force is barely possible, unless one has accurate optical 
data, extending over a wide range of frequencies around 
uj c , taken on the actual surfaces used in the experiment. 



III. A SUPERCONDUCTING CASIMIR 
EXPERIMENT 

In the light of the above considerations, we wondered 
if one can conceive a new type of experiment specifically 



devised to probe the magnitude of the contribution from 
the TE zero mode. As we see below, a superconducting 
Casimir apparatus is in principle a good tool to do that. 

Consider then two metallic plates made of a supercon- 
ducting metal, and imagine cooling them at a tempera- 
ture T\ < T c , where T c is the critical temperature. We 
suppose that the equilibrium temperature of the system 
is rapidly changed from the temperature T± to a higher 
temperature T% < T c . It is our aim to estimate the dif- 
ference AP between the Casimir pressures at the two 
temperatures: 



AP(a|T 2 , s;T u s) = P{a\T 2l s) - P(a|T 1)S ) 



(11) 



Let us explain the meaning of the new notation used in 
the above Equation. The notation is used to remind the 
two-fold dependence of the Casimir pressure on the tem- 
perature T. Indeed, by looking at Eq.JT]) we see that, 
on one side, P has an explicit temperature dependence, 
determined by the overall T-factor, in front of the sum 
symbol, and by the T-dependence of the Matsubara fre- 
quencies £/ = 2irkBTl/Ti. More importantly for our pur- 
poses, there is however also an implicit dependence of 
P on T, through the reflection coefficients, which may 
themselves depend on the temperature. The notation in 
Eq. (fTTjl stresses this fact, by making explicit reference to 
the dependence of P on the state of the plates, which can 
either be normal (n) or superconducting (s). We remark 
that probing the explicit T-dependence of the Casimir 
pressure, by letting pretty large changes of temperature 
(around 50 K), was the goal of the room-temperature ex- 
periment proposed in Ref. [10|. In that experiment, the 
temperature-dependence of the reflection coefficients was 
completely negligible. In our case the situation is much 
the opposite: the main expected cause of variation of 
P(a, T) is now the temperature dependence (in a certain 
frequency range) of the reflection coefficients of the plates 
in the superconducting state. 

After these remarks, we can go now to the main point: 
why should the proposed experiment be capable of telling 
us anything about the thermal Casimir problem in metal- 
lic systems? This is so, because the superconducting 
transition affects the optical properties of a metal only at 
frequencies ui smaller than a few times ksTc/h [20(. As a 
result, the magnitude of AP turns out to be almost com- 
pletely determined by the TE zero mode, and therefore 
it is very sensitive to the prescription adopted for it. We 
shall indeed show below that the two prescriptions con- 
sidered in the previous Section, lead to sharply different 
magnitudes for AP. 

Before we embark in the computation of AP, a key 
problem to address is to decide what prescription are we 
going to use for the reflection coefficient of the TE zero 
mode in the superconducting state. This is the subject 
of the next Section. 



IV. THE TE ZERO MODE IN A 
SUPERCONDUCTOR 

According to Lifshitz formula, which we recall is valid 
also for superconductors, in order to estimate the con- 
tribution of the TE zero mode to the Casimir pressure 
for a superconducting cavity, we need to determine what 
expression for the reflection coefficient ?te(0, kjj should 
be used in Eq.([T]), when a plate is superconducting. Ad- 
dressing this problem is the purpose of this Section. 

We have seen in Section II that there is a theoreti- 
cal uncertainty in the value of the reflection coefficient 

Tte(0, kj_) of a metallic plate, in the normal state, de- 
pending on whether we include or not, in the permittivity 
of the metal, the effect of ohmic dissipation. If dissipa- 
tion is included, we end up with the Drude prescription, 
Eq.([5]), if dissipation is ignored we instead obtain the 
plasma prediction, Eq.©. What about a superconduc- 
tor? For sure, static currents flow in a superconductor 
without any measurable dissipation, and so there should 
really be no room for ambiguity now: for what we know, 
in the dc limit, superconductors are strictly dissipation- 
less. This is reflected in the plasma-like form of the per- 
mittivity function e s (u>) that can be used to describe the 
response of superconductors to quasi-static electromag- 
netic fields, in the local London limit: 



,(«) = -[c/(A £ «)] s 



(12) 



In this Equation, Xl represents the penetration depth, 
that can be expressed in terms of the plasma frequency 
tip as Xl = (n/n s ) 1 ' 2 c/f2p, where n s /n represents the 
fractional density of superconducting electrons. The tem- 
perature dependence of n s /n is rather well described by 
the phenomenological law n s /n = 1 — (T/T c ) 4 [19[. If we 
insert e s (oj) into the Fresnel formula for the TE reflection 
coefficient at zero frequency we obatin: 
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for T <T C (13) 
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where we stressed that r^ is temperature dependent, 

because Xl is. The above form of r^(0,kx) is phys- 
ically plausible, and it correctly describes the Meissner 

effect. It should be noted that tvj?e(0, kx) has the same 
form as the plasma-model prescription r.j!g(0,kx)|pi for 
the normal state, apart from the replacement of H, P /c 2 
by 1/A|. In fact, taking account of the temperature de- 
pendence of Xl, r^g(0,kj_|T) provides a smooth inter- 
polation between rvj^j(0,kj_)| p i and 7Vj^j(0,kx)|Dr> as the 
temperature is increased from zero towards T c : 



r«(0,k x |T) 



rW( 0j k ± |T) 



41(0,kx)| P i for T/IWO, (14) 



r^(0,kx)| Dr = for T/T c - 1. 



(15) 



Now we are faced with the following question: depend- 
ing on what prescription we choose in the normal state, 
what do we do in the superconducting state? We can- 
not offer a certain answer to this question, but the fol- 
lowing choices appear to us as reasonable. Consider first 
the Drude prescription. Since the superconducting phase 
transition is of second order, the reflection coefficients 
must be smooth across T c . Now, the Drude prescription 
asserts that tv-j^O, kj_) is zero, and therefore we should 
fix the prescription in the superconducting phase in such 



», 



Tc. 



a way that r^(0,]f.±\T) approaches zero for T 

Of course, we could achieve this by taking rvjfj^O, kj_|T) 
identically zero for T < T c . However, such an ansatz 
looks to us physically wrong, because it is in conflict with 
the Meissner effect. The prescription that we choose then 
is to use Eq. (|f3[) . with Xl temperature dependent, as our 
Drude recipe in the superconducting phase: 



0) 

r TE 



(0,kx|T)|i 



0) 

r TE 



(0,kjT) forT<T c . (16) 



This choice ensures smoothness of the reflection coeffi- 
cient at the phase transition, while describing correctly 
the Meissner effect. Importantly, within this prescrip- 
tion, the reflection coefficient is temperature dependent 
in the superconducting phase. 

Consider now the plasma prescription. Again, the re- 
flection coefficient must be smooth through the transi- 
tion. If we further make the plausible assumptions that 
the reflection coefficient cannot decrease as we reduce the 
temperature, and that the plasma frequency sets up an 
upper limit on the quantity c/Xl, we are led to conclude 
that, within the plasma prescription, the reflection coef- 
ficient ripg(0,kj_|T) has to be temperature independent 
and equal to its normal-state value: 



r TE 



(0,k x |T)| p i = r^(0,kj.)| pl forT<T c . (17) 



It should be noted that, with this prescription, the reflec- 
tion coefficient does not change at all across the transi- 
tion, and it is temperature independent in the supercon- 
ducting phase. This feature marks a crucial difference 
between the plasma prescription and our Drude prescrip- 
tion in Eci. tfTS]) . 

Having defined our prescriptions for the TE reflection 
coefficient in the superconducting state, we are ready now 
to proceed with the computation of AP according to the 
two approaches. 



V. COMPUTING AP 

In this Section, we undertake the computation of AP. 
We shall consider two possible geometries: the plane- 



J 



parallel configuration, and the sphere-plate configura- 
tion, which is the one usually adopted in current ex- 
periments. For either geometry, it is useful to consider 
two different setups for the system: the first one consists 
of two superconducting plates, while the second setup 
has just one superconducting plate, the other plate be- 
ing made of a normal metal. To be definite, we shall 
consider Nb as our superconducting material, and Au as 
the normal metal. The choice of Nb is due to the fact 
that, among the classic metallic superconductors, it is 
the one with the highest critical temperature T c = 9.2 
K, a feature that will be seen to ensure the largest differ- 
ence between the Drude and the plasma predictions for 
AP. It is assumed that the Nb plates have a thickness 
much larger than the penetration depth Xl of magnetic 
fields. In this Section we shall only consider the plane- 
parallel case. The sphere-plate case will be treated in 
Sec. VI below. 

It is convenient to split Eq. (fTTjl in a way analogous to 
Eq. ©: 



AP = AP 



(TE) 



AP, 



(TM) 



APi, 



(18) 



where the meaning of the symbols is obvious. Leaving 
aside for the moment APq , we consider first AP 
and APi. Since ?"TM(0,k_L) is one for metals, no matter 
whether normal of superconducting, we surely have: 

Ap(™\a\T 2 ,s;T 1 ,s) = Ap(™\a\T 2 ,n;T 1 ,n) . (19) 

Consider now APi. Recall that APi includes only the 
contributions from non vanishing Matsubara modes. The 
key experimental fact to consider now is that the opti- 
cal properties of superconductors are appreciably differ- 
ent from those of the normal state (right above T c ) only 
for frequencies less than a few times ksTc/h [20]. Since 
£i(T 2 ) ~ £i(Ti) ~ £i(T c ) = 2nk B T c /h, we see that al- 
ready the first non- vanishing Matsubara mode is over six 
times larger than kBT c /h. Therefore, in the range of 
frequencies that is relevant for APi, the superconductor 
can be described, at all temperatures below T c , by the 
same set of reflection coefficients of the normal state, at 
temperatures slightly above T c (this consideration applies 
even more to Au, of course). This implies at once that 



APi (a \T 2 , s; T x ,s) = APi(o \T 2 , n; T u n) 



(20) 



At this point, it is convenient to add and subtract from 
Eq.flTU) the quantity Ap[ ) TE) (a \T 2 , n; T x ,n). In this way, 
we get 



AP(a|T 2 , S ;T 1 ,s) 



[AP( TE) (a |T 2 , S ;Ti,s) 



AP<, TE) (a|T 2) n;Ti,n)] 



[AP^\a\T a ,n;T u n) + AP^(a\T2,n;T 1 ,n) + M> 1 (a\T 2 ,n;T 1 ,n)]. 
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(21) 



Now we note that the sum of three terms on the second 
line of the above Equation is nothing but the quantity 
AP(a \T 2 , n; T\, n) that represents the change of Casinhr 
pressure when the temperature is changed from T\ to 
T 2 . in a normal metal system, with temperature inde- 



J 



pendent reflection coefficients. This quantity has been 
computed already in Ref . [lfjl| . both for the Drude and 
the plasma prescriptions, and we shall later use the for- 
mulae derived in that paper to estimate it. Substituting 
AP(a \T 2 , n; T\,n), we may then rewrite Eq.([2lTl as: 



AP(a\T 2 ,s;T 1 ,s) = [AP^ E) (a\T 2 ,s;T l , S )-AP^ E) (a\T 2 ,n;T l ,n)}+AP(a\T 2 ,n;T 1 ,n) 



r 



(22) 



The above Equation represents the first key result of this 
paper, because it shows, as announced earlier, that super- 
conductivity has an effect on the Casimir pressure only 
via the TE zero mode. At this point we consider sepa- 
rately the plasma and the Drude prescriptions. 



A. The plasma prescription 

According to our plasma prescription for the supercon- 
ducting state, Eq. lfTT)) . the reflection coefficient for the 
TE zero mode in the superconducting state has precisely 
the same expression as in the normal state. This being 
so, we obviously have: 



,(TE) 



(TE), 



AP^^(a|T 2l s;T 1 ,s)| pl -AP^ J1J (a|r 2 ,n;T 1 ,n)|pi = 0, 

(23) 
both for the Nb-Nb and the Nb-Au setups. When this 
formula is used in Eq. (|2"2")l . we see that, in either setup, 
superconductivity has no effect on AP: 

AP{a \T 2 , s; T 1: s)| pl = AP(a \T 2 , n; T u n)| pl . (24) 

The quantity AP(a \T 2 ,n;Ti, n)\ p \ has been computed 
in Ref. [10], for the case of two plates made of the same 
metal. This case is good enough for us, also in the case 
of the Nb-Au setting, because of the small difference be- 
tween the plasma frequencies of Nb (8.7 eV/h) and Au (9 
eV/h). In our computations, from now on, we shall than 
take for both metals the common value flp = 9 eV/h. 
Then, our AP(a \T 2 , n; T\, n)| p i coincides with minus the 
quantity AF pp of Ref.[l0J (the minus sign is due to the 
fact that we consider positive forces to correspond to an 
attraction, while Ref.[T(| uses the opposite convention). 
Therefore, we can rewrite Eq. (f2"4"|) as 



(25) 



AP(a\T 2 ,s;T 1 ,s)\ pl = -AF p: 



We can obtain an estimate of the effect by using the 
following perturbative expression for AF pp of Ref. [1C|| . 
that holds for low temperatures kBTa/(hc) <C 1, and for 



not too small separations cj c /f2 p <C 1: 
AF w (o,Ti,T 2 ) = -A^F pp (T x ,T 2 )A^F pp (a,T u T 2 ), 
where 



A(1)P frp T\- ff ^b( T 2 T l) 

A^ Fpp (T u T 2 )- 45? .3 c3 



(26) 

(27) 



and 



T{T 2 



I? 



•n 



7r d a Ti + T 2 \ 

where ksT c ff = he/ '(2a) and 6 — c/Cl p . The above formu- 
lae show that AP(Nb - Nb)| p i and AP(Nb - Au)| p i are 
both extremely small. For example, upon taking Pi = 5 
K, T 2 ~ T c , and a = 100 nm, we obtain: 



AP 



pi 



1.4 x 10~ 9 Pa. 



(29) 



The conclusion that we draw is that, with the plasma 
prescription, the change with temperature of the Casimir 
pressure is unmeasurably small, both in the Nb-Nb and 
in the Nb-Au setups. 



B. The Drude prescription 

We now consider the Drude prescription. First, we 
analyze the Nb-Au setup. Since Au is always a normal 
metal, its Drude reflection coefficient for the TE zero 
mode is zero, at both temperatures. Therefore, since 
Lifshitz formula implies only products of reflection coef- 
ficients for the two plates, it follows that the TE zero 
mode gives no contribution, in the Nb-Au setup: 



AP< TE) (a|r 2)S ;Ti, S )|g 
AP< TE) (o|T 2> n;r 1> n)|g 



Nb-Au 



Nb-Au 



0. 



(30) 



Upon substituting these formulae into Eq. (|22|) we obtain 



AP(a\T 2 ,s;T 1 ,s)\^ 



Nb-Au 



AP(a\T 2 ,n;T 1 ,n)\^ 1 . 



Nb-Au 



(31) 



Again, we find that superconductivity has no effect on 
the change of Casimir pressure. The Drude value of the 
quantity AP(a |T2,n;Ti,n)| Dr _ u can also be found in 
Ref. [KJ], from which we obtain the estimate: 



Nb-Au 



C(3)*b(T 2 -Ti) 



= -AF, 



i>i> 



&ira 3 



1-6- 



24- 



(32) 



The expression on the right hand side of the above 
Equation, differs from the analogous expression obtained 
within the plasma prescription, Eq. (|25p by the presence 
of the extra term on the second line. At our cryogenic 
temperatures, the absolute value of the latter term is 
larger than |AP PP | by several orders of magnitude, for 
changes of temperature of a few degrees. We therefore 
see that heating the system, the Casimir pressure de- 
creases by an amount proportional to the the temper- 
ature change. We remark once again that this change 
of Casimir pressure is not an effect of superconductiv- 
ity, and it only arises from the explicit T dependence of 
the Lifshitz formula. Indeed, an analogous phenomenon 



was found to occur in the room temperature setting of 
Ref. [HI. 

We finally consider the Nb-Nb setup, the most inter- 
esting case indeed. Using the reflection coefficient in 
Eq. (fi~6|) . we easily find 



,(TE) 



Nb-Nb 



AP^(a|T a ,«;T 1( 5)|£ 



P (TE) (a, T 2 ; c/\ L ) - P (TE) (a, Ti; c/X L 



>(TE) 



Ajy^ColTa.njTi.TOIS 



Nb-Nb 



0. 



(33) 



(34) 



where, in the first of the two above Equations, 
P (a, T; c/Xl) denotes the magnitude of the TE zero 
mode, that results after we plug into the I = term of 
Lifshitz formula the expression the Drude ansatz for the 
TE reflection coefficient for the superconductor, Eq. (fl"6|) . 
As for the quantity, AP(a \T2,n;Ti,n)\x>r, its magnitude 
is given by the same expression on the right hand side of 
Eq. ([32]) . Collecting everything together, Eq. (|22|) gives: 



J 



AP(a|T 2 , S ;T 1 , S )|S r c 



Nb-Nb 



[P (TE) (a,T 2 ;c/A L ) 



>(TE) 



(a,Ti;c/A L )]-AF, 



pp 



C(3)fcB(T 2 -Ti) 

87T0 3 



6* 



S 2 
a a z 



J351 



In Fig.l we show plots of AP(a \T 2 , s; T x , s)|^" Nb (solid 
line) and AP(a \T 2 , s; Ti, s)|gJf~ Au (dashed line), both ex- 
pressed in mPa, for a = 150 nm, as a function of the 
temperature Ti (in K), for T 2 ~ T c . In Fig. 2 we show 
the plots of the same quantities, this time as functions 
of plates separation a (in microns), for T\ = 5 K and 
T 2 ~ T c . Note that, contrary to what we found with 
the plasma prescription, our Drude prescription predicts 
negative changes of pressure, i.e. the Casimir pressure 
decreases as one goes from the superconducting state to- 
wards the normal state. Moreover, it should also be noted 
that the change of pressure is larger in the Nb-Nb setup, 
than in the Nb-Au setup. However, the really striking 
finding is that, for both setups, the Drude approach pre- 
dicts changes of Casimir pressure, whose magnitudes are 
around five or six orders of magnitude larger than the cor- 
responding changes predicted by the plasma approach! 
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FIG. 1: Plots of AP(Nb - Nb)| D r (solid line) and AP(Nb - 
Au)|Dr (dashed line) (in mPa) for two parallel plates at fixed 
separation a = 150 nm, versus temperature Ti (in K), for 
fixed T2 — T c . 



VI. THE SPHERE-PLATE CASE 



Now we consider the configuration of a sphere placed 
above a plate, which is the one used in most of the current 
experiments. We let R the radius of the sphere and a the 
minimum sphere-plate separation. As it is well known, 
for small separations, i.e. for a/R <C 1, the Casimir 
force F sp between a sphere and a plate can be obtained 



by means of the so called proximity force approximation 
(PEA) [2J. The relative error introduced by this approx- 
imation is of order a/R [2], and therefore it is only a 
fraction of a percent in typical experimental situations. 
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FIG. 2: Plots of AP(Nb - Nb)| Dr (solid line) and AP(Nb - 
Au)|Dr (dashed line) (in mPa) for two parallel plates versus 
separation a (in microns), for 2\ = 5 K and T2 ~ T c . 



where a plus sign for the force again corresponds to 
attraction. Analogously to what we did in the plane- 
parallel case, we consider two distinct setups: in the first 
one, both the plate and the sphere are made of Nb, while 
in the second setup the sphere is made of Nb and the 
plate of Au (or viceversa). For either setup, we consider 
then the change in the sphere-plate force 



AF ps (a|T 2 , s; T x , s) = F ps (a|T 2 , s) - F ps (a\Ti, s) , (37) 



The PFA expression for F ps is: 



k B TR ^ , 



2~ 



IX) 



d 2 k. 



q=TE,TM 



log(l 



_ r (l) r (2) e -2a ?! 



(36) 



that occurs when the equilibrium temperature of the 
system is increased from from T\ to T 2 (the notation 
adopted here is analogous to that used in the plane par- 
allel case). By repeating the reasonings that led us to 
write Eq. (|22p , we can obtain the following expression for 
AF ps (a\T 2 ,s;T 1: s) 



AF ps (a \T 2 ,s;T 1 ,s) = [AF < ^l (a \T 2 ,8 ] T u a)-AF™(a\T 2 ,niT 1 ,n)] + AF„(a\T a ,n;T 1 ,n), 



(38) 



r 



where again the meaning of the symbols is obvious. By 
similar steps that led us to write Eq. ([24|) . one can show 
that, both for the Nb-Nb and the Nb-Au setups, the 
plasma prediction for the change of force in the sphere- 
plate geometry is: 

AF ps (a|T 2! s;T 1 ,s)| pl = A^ ps (a|T 2) n;Ti,n)|pi. (39) 

Once the small difference between the plasma fre- 
quencies of Nb and Au is neglected, we can use the 
explicit expressions provided in Ref.[10j to evaluate 
AF ps (a |T 2 , n; T\, n)| p i (recall that our sign convention 
for the Casimir force is opposite to that of [10|): 



AF ps (a|T 2 ,n;Ti,n)| p i = 
RA^F ps {T x ,T 2 )A^F ps {a,T Xi T 2 ), 



(40) 



where 



A^F ps (T 1 ,T 2 ) = 



C(3)*|^ ^w, 



2 2 



h c 



{T 2 -T{){Tt + Ti), (41) 



and 



A^F p8 {a,T x ,T 2 ) = (1 



TiT 2 



n 



n 



Ti+T 2 



45C(3) Teff 



,2 S - 
a 

4* 
a 



(42) 



Analogously to what we found in the plane-parallel case, 
the effect predicted by the plasma prescription is un- 
measurably small. For example, for a sphere of radius 
R = 200 /im, at a distance a — 150 nm from a plane, for 
T\ = 5 K and T 2 ~ T c the above formulae give a change 
of Casimir force of about 5.3 x 10~ 19 N. 

We consider now the Drude prescription. Repeating 
the steps that led to Eq. (|4"3"l) , we obtain for the Nb-Au 
setup the Equation: 

AF pB (o|T 2 ,*;Ti, a )|^- Au = AF pil (o|T 2 ,n;T lj n)|S- Au . 

(43) 
An explicit formula for AF ps (a \T 2 , n; T 1: n)\^- Au is 
given in Eq.(10) of Ref.0: 

AF ps (a|r 2 ,n;T 1 ,n)|^- A 

i?fcsC(3) 



8a 2 



(T 2 - Tx) 



RA^F ps A^F ps 

S S 2 \, N 

1-4- + 12-J 1(44) 



The first term on the right hand side of the above Equa- 
tion is completely negligible with respect to the second 
term. We therefore see that the Drude model predicts 
a decrease of Casimir force proportional to T 2 — T\. We 
repeat that this effect has nothing to do with supercon- 
ductivity, and it is again analogous to what was found in 
Ref . ( [lOJI ) (see comments following Eq. ([3^1) ). 

It remains to discuss the Nb-Nb setup. By proceeding 
in the same way as we did to derive Eq. (|35p . one can 



-2 



rt 



Em 




prove the following Equation 



FIG. 3: Plots of AF sp (Nb - Nb)| Dr /-R (solid line) and 
AF sp (Nb - Au)| Dr /-R (dashed line), in units of 10" 10 N/m, 
versus Ti (in K), for T% — T c and for fixed sphere-plate sepa- 
ration a — 150 nm. 
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FIG 

AF sp (Nb - Au)| Dr /-R (dashed line), in units of 10" 10 N/m 
versus sphere-plate separation a (in microns), for T\ — 5 K 
and T 2 ~ T c . 



AFp 8 (o|r 2 ,a;T 1)a )|^ 



Nb-Nb 



[F^(a,T 2] c/X L ) 



F^(a,T i; c/X L )} 



RA^F ps A^F ps 



^® (T, - TO f 1-4^ + 124 

8a z V a or 



(45) 



r 



where Fi ps J (a, T; c/ A_l) denotes the TE zero mode con- 
tribution to the sphere-plate force (the 1=0 term for TE 
polarization in Eq. ([36| ). as it results after we use our 
Drude prescription for the reflection coefficient of the 
superconductor, Eq. tfTo) . In Fig. 3, we plot AF ps (Nb — 
Nb)| Dr /i? (solid line), and A.F ps (Nb- Au)| Dr /i? (dashed 
line), in units of 10~ 10 N/m, as a function of T\ (in K), 
for T 2 ~ T c and for a plate-sphere separation a = 150 
nm. In Fig. 4 the same quantities are plotted versus plate- 
sphere separation a (in microns), for fixed temperatures 



Ti = 5 K and T 2 ~ T c . Note that the changes of force 
are always negative, indicating that the Casimir force 
decreases as the temperature of the superconducting sys- 
tem is increased from T\ to T c . As we see from Fig. 3, 
for a sphere-plate separation of 150 nm and for a sphere 
of radius R = 200 fim, the Drude model predicts for the 
Nb-Nb setup a change in the force around —0.8 x 10~ 13 
N, which is over five orders of magnitude larger than the 
plasma model prediction! 



10 



VII. CONCLUSIONS AND DISCUSSION 

One of the most important unresolved theoretical 
problems in the theory of dispersion forces is that of de- 
termining the thermal correction to the Casimir force 
between two conductors. The source of the difficulties 
stems from the ambiguity in the correct value for the 
TE zero-mode contribution to the force, corresponding 
to quasi-static magnetic fields. In this paper, we have 
proposed a new Casimir experiment with superconduct- 
ing cavities, specifically devised to probe the TE zero- 
mode. The proposal involves measuring the change of 
Casimir force that occurs in a cavity with one or two su- 
perconducting plates, as the temperature of the device 
is increased towards the normal state. The interest of 
performing difference force measurements to probe the 
features of the thermal correction to the Casimir force 
has already been stressed in Ref.[l(|, which considered 
measuring the difference in the Casimir force at two dif- 
ferent temperatures in cavities made of ordinary metals 
at room temperature. Our proposal is very much in the 
same spirit as Ref. 10], with the significant difference that 
the change in the Casimir force determined by a small 
change of temperature is much enhanced (in the Drudc 
case) by superconductivity of the plates. The need of 
smaller temperature changes, than those required in [lfj, 
might make it easier to achieve them in an experiment. 

It is perhaps the case to comment further on the pos- 
sible advantages and drawbacks that are implicit in the 
proposed scheme. A potential advantage is that, being 
a difference force measurement, it might be possible to 
achieve a better sensitivity than in an absolute force mea- 
surement. This of course holds true only provided that 
systematic errors do not change too much when the tem- 
perature of the apparatus is varied. Importantly enough, 
this is the case with regards to the theoretical uncertainty 
arising from insufficient knowledge of the optical data of 
the plates, that was estimated in Ref. [18| to be easily 
as large as five percent in absolute force measurements. 
This uncertainty is completely irrelevant in our scheme, 
because as we showed in Sec. Ill, the magnitude of the 
effect is determined solely by the TE zero mode. More 
delicate is the case of the systematic errors arising from 
residual electrostatic attractions between the plates [17| , 
which constitute a concern in all Casimir experiments. If 
their magnitude is temperature dependent in the super- 
conducting state, it might be necessary to perform elec- 
trostatic calibrations of the apparatus at the considered 
temperatures. 

One of the main experimental issues to address is to 
find a good mean of varying the temperature of the 
plates. Perhaps, a convenient way to do this is by il- 



luminating the plates with short laser pulses, as sug- 
gested in Ref. flOJ] . Very likely, this method of heating 
the plates will reduce to a minimum unwanted thermal 
expansions/contractions of parts of the apparatus that 
would inevitably result, were we to warm up the whole 
system instead. Such expansions would clearly easily al- 
ter the separations of the plates by a few nanometers, 
thus rendering much more difficult a comparison of the 
Casimir force between the two temperatures. 

Finally, we would like to comment on the prospects 
of the experiment described in this paper of being ac- 
tually capable of discriminating between the contradic- 
tory theoretical approaches to the description of the ther- 
mal Casimir force between real materials. According 
to our computations, the two approaches predict strik- 
ingly different magnitudes for the temperatures changes 
of Casimir force in the superconducting state, the plasma 
approach predictions being always five or six orders of 
magnitude smaller than those of the Drude approach. 
While the changes of force resulting from the plasma ap- 
proach are unmeasurably small, the effect predicted by 
the Drude approach is in fact large enough for an ex- 
perimental detection to be hopeful. For example, in the 
sphere-plate geometry, which is the most widely used ge- 
ometry in current experiments, if both the sphere and the 
plate are made of Nb, for a sphere radius of 200 /im, and 
a sphere-plate separation of 150 nm the Drude approach 
predicts a drop in the Casimir force of 0.8 x 10~ 13 N, 
as the system is heated by two-three degrees. Force os- 
cillations slightly larger than this one have been demon- 
strated to be measurable by means of an atomic force 
microscope, with an error about 20 percent at 95 per- 
cent confidence level, in a recent measurement of the so- 
called lateral Casimir force [2l[ . If the same accuracy can 
be achieved at cryogenic temperatures, there are good 
chances that the experiment proposed here might actu- 
ally be able to discriminate the alternative approaches to 
the description of the thermal Casimir force between real 
metals. 
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